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Using one loop functional RG we study two problems of pinned elastic systems away from their 
equilibrium or steady states. The critical regime of the depinning transition is investigated starting 
from a flat initial condition. It exhibits non trivial two-time dynamical regimes with exponents 
and scaling functions obtained in a dimensional expansion. The aging and equilibrium dynamics of 
the super-rough glass phase of the random Sine-Gordon model at low temperature is found to be 
characterized by a single dynamical exponent z ~ c/T, where c compares well with recent numerical 
work. This agrees with the thermal boundary layer picture of pinned systems. 
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Disordered elastic systems offer many experimental re- 
alizations and are also of theoretical interest as prototype 
models for glasses induced by quenched disorder. The 
competition between the structural order and substrate 
impurities results in pinning, complex ground states, bar- 
riers and ultra slow glassy dynamics. While ground state, 
equilibrium dynamics, and driven steady state properties 
have been much studied theoretically, much less is known 
about the dynamics before the steady state is reached, 
or about the aging dynamics. If universality is shown 
there, it would be of high interest for numerous exper- 
imental systems, e.g. magnetic domain wall relaxation 
0, superconductors contact line depinning |2iJ, den- 
sity waves 0. 

Numerical studies of glassy dynamics are hampered by 
high barriers in configuration space resulting in ultra-long 
time scales making comparison with theory uncertain. In 
some cases however faster, but still interesting dynamics 
occurs. One is zero temperature T = driven dynam- 
ics near the depinning transition where barriers disap- 
pear 0. Recent theoretical progress has been achieved 
there. Functional renormalization group (FRG) studies 
give more precise and consistent predictions , corrobo- 
rated by powerful new algorithms which allow for excel- 
lent determination of steady state exponents 0, 0] . One 
aim of this paper is to extend the FRG to dynamics away 
from steady states, and to show that interesting univer- 
sal two-time dynamics (analogous to the aging dynamics 
in non driven situations) also occurs near depinning. We 
predict new exponents and scaling functions in the criti- 
cal regime. 

A second case where T > dynamics has been inves- 
tigated is the "marginal glass" phase exhibited by topo- 
logically ordered 2D periodic systems as captured by the 
Cardy Ostlund (CO) model Barriers there grow only 
logarithmically with size allowing for precise numerics 
[lOj. The equilibrium and aging [1^ dynamics were 
studied using Coulomb gas RG methods, but only near 
the glass transition temperature Tg A numerical 

study has confirmed some of the RG predictions, and in 



addition has explored the full temperature regime jl4l |. 
The dynamical exponent was found to diverge as z 1 /T 
at low T. One aim of this paper is to show this result 
within a simple one loop FRG, and to obtain detailed 
predictions for the aging regimes at low T . This study, 
together with the rather good agreement with numerics, 
is also important as indirect evidence for a recent hypoth- 
esis that a "thermal boundary layer" (TBL) in the field 
theory controls the activated dynamics of (more strongly) 
pinned manifolds (isf. 
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FIG. 1: (a) solid line: initial flat configuration, dashed and 
dotted lines: configurations of the string at time t' and t re- 
spectively, (b) Domain of validity, in logarithmic scale, of our 
approach near the depinning transition. The grey circle rep- 
resents the small time region sensitive to microscopic details. 

The overdamped dynamics of a single component elas- 
tic manifold of internal dimension d, parametrized by a 
height field u{x,t) = Uxt is described by 

r]9tUx,t = c\''^Ux^t + Fix,Ux,t + vt) + f -r]v + ^{x,t) (1) 

where = 0, (C(x, t)^(x', t')) = 2Tr]S'^{x - x')S{t - 
t') is the correlator of the thermal gaussian noise, 
F{x,u)F{x' ,u') = A{u — u')S'^{x — x') the second cu- 
mulant of the quenched random pinning force and v is 
the average velocity. We denote the small scale cutoff 
size a A^^. In this paper we consider a flat initial con- 
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figuration Uxt=o = 0. Denoting Uqt the spatial Fourier 
transform of Ux^t, we focus on the correlation C^^, and 
the connected (w.r.t. the thermal fluctuations) correla- 
tion C«,: 



Cj'j, = {UqtU-qt') , Cj^t' = {UqtU-qt') - {Uqt){u^qt') 



and the response 7^?^, to a small external field 



qt' 



5{Uqt) /Sh^qf 



t > t' 



(2) 
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where — and (..) denote respectively averages w.r.t. dis- 
order and thermal fluctuations. In the following i?'^, = 
7y~^e~' (*-*)/'; (respectively C^^,) denote the bare re- 
sponse function (respectively correlation), i.e. in the ab- 
sence of disorder. We now specify the time regime t, t' . 

In the first situation studied below, i.e. an interface 
very near the depinning threshold, we set T = in 
Numerical simulations of a manifold of internal size L are 
typically performed on a cylinder, i.e. with periodicity 
in u ^ u + W with W ^ L''. After a time empirically of 
order fixed number of turns, i.e. according to 

scaling (z the dynamical exponent), the system reaches a 
unique global time-periodic steady state which has been 
extensively studied (3, IS] • Here we are interested instead 
in the dynamical regime before this steady state is estab- 
lished, i.e. in the limit of times t, t' very large compared 
to microscopic time scales a^, and such that: 



t, t' <^L\e 



(4) 



where ^ ~ \ f — fc\~'^ is the scale above which interface 
motion becomes uncorrelated, i.e. we mainly focus on 
the critical regime of highly correlated avalanche motion 
(Fig. 1). In the computation of the correlations Q and 
response we will be interested in the scaling limit 
g/Ao ^ 1, keeping the scaling variables 



w = q\t-t') 



t/t' 



(5) 



fixed. In the context of the depining transition, the do- 
main of validity of our approach is depicted in Fig. 1. 

To study the second situation, i.e. the relaxational 
dynamics of the c? = 2 Cardy-Ostlund model, we set / = 
V — in (Q, A(u) being a periodic function of period 
unity. The super-rough glass phase ^ll'l for T < Tg is, 
within RG, described by a line of finite temperature fixed 
points (FP). Recently, we have obtained analytically ^3 
Cf^i © under the scaling form: 



T 

C, - — 



(6) 



where Fq(w, u) is a universal |l8| scaling function, which 
was recently confirmed by numerics |l4l | in a wider range 
of temperature. And although the numerically measured 
exponents zco and A were found to be in good agree- 
ment with one loop RG predictions near Tg, significant 



deviations were found to occur at lower T. In addition, 
the (connected) structure factor is obtained from Eq. 
© in the limit w — > 0, w — > 1 keeping ■w/{u—l) = q^-'^°t 
fixed. Therefore the dynamical exponent zqq associated 
to equilibrium fluctuations coincides with the one asso- 
ciated to non- equilibrium relaxation, which was actually 
numerically computed in Ref . |l4l | . 

In both cases Eq. (Q) is studied using the standard dy- 
namical (disorder averaged) MSR action. Correlations 
and response Q are obtained as functional deriva- 
tives of the dynamical effective action r[M, u]. It is per- 
turbatively computed using the Exact RG equation 
(ERG) associated to the multi-local operators expansion 
introduced i n and extended to non-equilibrium 

dynamics in The ERG equations are obtained by 

varying an infrared (large scale) cutoff A; = Aqc"' in- 
troduced in the bare response and correlation functions 
{R ^ Ri, C ^ Ci). The information about non equilib- 
rium dynamics is contained in the interacting part of F: 



int = I iUxtFlt [u^ 



iuxtiuxt'^itt'{u^t - Uxt') 
(7) 

where only the solution of the ERG equation to lowest 
order in A = A; and to one loop is needed here. It reads: 



SUxt' 



Rf-,°A"{uxt-Uxt') 



(8) 



AuAu) = A(^) -t- A"(u)(cr,^" - icr,r° - \cr,r°) 

As in we need however the FRG equation for the 
"statics" part to one loop and next order. It has the 
standard form: 



diA{u) = (e - 2C)A(u) + CuA'(u) + fA"{u) 
-l[(A(u)-A(0))Y' 



(9) 



where one defined the rescaled, dimensionless disorder 
Ai{u) = SdA^^Aiiu) and temperature fi = SdAf~'^T/c 
with e — 4— d. The response TZ^^, and correlation function 
C^j, at the fixed point {I 00) are given, to one loop, by 
the equations: 



f dh^tt,n%, - f dt,i:u,nl,, (lo) 

Jo Jo 



Jti Jti Jti 

SF u ^^^^ 

where the fixed point self-energy is Stf = ^^'^^^"^ \u=o 
and the disorder-noise kernel Du' — Aiu'{u — 0) are 
local in space to this order of computation. 

We first apply the above equations to the case of the 
depinning transition, just above threshold / = f^. In 
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that case one further introduces a small finite velocity 



in the above equations (shifting u^t 



Uxt + vt 



evry where). Eq 0, setting T = reaches the standard 
one loop FP for depinning transition with ^^ep = e/3- 
This FP function being non-analytic at m = 0, this results 
in Y,tv = i?ftp° A"(0+) using the limit v = 0+. We can 
now solve the equation for 7^'^, (|10() . perturbatively in 
the disorder, as in |l2l |. and obtain a solution consistent 
with the scaling form 



(12) 



with z — 2 = — A*"(0+) and the novel exponent Or asso- 
ciated to large time off equilibrium relaxation: 



9n = ^^A*"{0+) = -'- 



(13) 



Fii{w,u) = Fii{w) is a universal [l3 | scaling function, 
whose expression is given at one loop order by: 

Fii{w) = e-"' + ^y^((w-l)Ei(ii;)e-™ + e-'^-l) (14) 

where Ei (w) is the exponential integral function, with the 
large w power law behavior Fji{w) oc w~'^ . This one loop 
scaling form q^~'^Ffi{q^{t — t')) can be written as the 
Fourier transform, w.r.t time variable of 1/(9^ -I- E(iw)) 
wfth Yj{s) = s^/z + s _|_ 0{s^). Such scaling forms lf]l|l 
arise in the context of critical points Solving Eq. 

(Il()|l for any finite Fourier mode q, one obtains the local 
response function 7?.f(7° 



(t_t')l + (d-2)/z 



(15) 



where A^,A\^ are non universal, Ao-dependent, ampli- 
tudes (the logarithmic corrections coming from the large 
w behavior of Ffi{w,u)). At this order, 11511 is compati- 
ble with local scale invariance arguments [20j . Note how- 
ever that at this order (|15|l could also be written as j2l] | 
n^=o ^ An {t/t'f" {t - <')-(!+«) with a ^ l + (d-2)/z : 
clarifying this point requires higher order calculations, 
left for future investigations. 

The correlation function is obtained by solving per- 
turbatively the equation Hll|l . at T = (thus the bare 
correlation C'f^^ and the connected one C'j, vanish). One 
finds in perturbation to one loop and lowest order that 
it is consistent with the scaling form: 



-(<i+2Cdep) 



Oc 



Fc{q^t~t'),t/t') (16) 



Fciw,u) = A*(0)(1 - e"'"'"/("^i))(l - e-"'/("-i)) 

where a priori 9c = 0{e) which, as Fc{w,u) is already 
of order 0{e), requires a second order calculation [2^ . 



We now focus on the relaxational dynamics of ran- 
dom periodic system, at low temperature. We first start 
by deriving the low T expression of the equilibrium dy- 
namical exponent zco- It is given by the one loop FRG 
equation © together with di Inrji = — A"(0) = zco — 2. 
We will specialize to d = 2 (we remind ^2 ~ l/(27r)) thus 
e ~ 2 and therefore Ti = T = T/ (27rc) is not flowing. In 
that case, there is a line of fixed points Ay(it) indexed 
by temperature, analyzed in Ref. |23| . for T <Tg. The 
transition ai T = Tg can be analyzed in the Fourier rep- 
resentation A(u) — A„ cos(27rnM), the linear part 
of Eq. © being 9/A„ = 2(1 - ^n2)A„ + ©(A^), where 

Tg = c/tt (i.e. 2 - (27r)2fg = 0). The transition cor- 
responds to the lowest harmonic becoming relevant. At 
any temperature the fixed point solution of ^ can be 
written in that case: 



u = (3/e)G(T',f -f A*(0) - A*(m)) 



(17) 



4y± = 3A*(0) -I- T ± V3(3A*(0) + r)(A*(0) + 3r) 



where G{a,b) = ydy/^Jiy - T){y - y^){y - y+), and 
y- < Q <f < y+. The condition i = G'(f,y+) yields 
A*(0) as a function of T. As T ^ this solution 
converges to the zero temperature solution Ay^g(u) = 
At=o(0) - §"(1 - u), and A^^o(O) = e/36. From the 
FRG equation at u = one has the exact relation for 
all T: 



f A*"(0) = -eA*(0) 



(18) 



It implies that as T ^ 0, A*"(0) eA^^o(0)/f which 

gives, using T = T/{2Tr^Tg) the one loop estimate: 



Zco 



2n^ T, 
~9 



a 
T 



2.19^ 



(19) 



We now compute the response T?,'^, (PJ and the corre- 
lation C^(, (O for the relaxational dynamics defined by 
Eq. By solving to first order in the disorder the 

equations I|1UI11|) for the present case, one obtains, in 
the limit g/Ao ^ 1 keeping q^{t — t'), t/t' fixed, the same 
result as for the the depinning H121 1161) with the substi- 
tution of the exponents z and by zqq given by Eq. 
(|19|l and Ocq given to one loop, as T — > by 



9 T T 



(20) 



For the present case where T ^ 0, the connected corre- 
lation function C^^, is non zero. It is given by an equation 
exactly similar to Eq. Hll|) with the substitution of Dte 
by Djj, given to one loop by D^^, = A"(0)Cj^^°. By solv- 
ing perturbatively the equation for Cj' / , we find a solution 
consistent with the scaling form given in Eq. ((HJ where 
Zco is given by lfTO|l and A = d to the order of our cal- 
culation, in good agreement with the numerics |l4l | and 
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FIG. 2: l/zco as a function of T /Tg. The square symbols are 
the resuhs of the numerical simulation of Ref . jl^ . The curves 
FRGl and FRG2 correspond to our one loop FRG estimate 
as explained in the text. For comparison, we have plotted in 
dashed line the result of the one loop RG calculation in the 
vincinity of Tg llj. 



independently of w, a consequence of (|^ . Notice also 
the identity, given the value of zco ifT^ - obtained here 
up to one loop, zco = This relation was also found 
in the vicinity of Tg and is consitent with numerical 
simulations at low T |l4| . 

In conclusion we have defined and computed new uni- 
versal exponents and scaling form for driven interfaces 
near the depinning transition. We also showed that the 
one loop truncation of the FRG yields a good approxima- 
tion to numerics for low T aging dynamics of the pinned 
periodic manifold in c? = 2. Further numerics near depin- 
ning and investigations of other predictions of the TBL 
picture (e.g. barrier fluctuations), together with more 
precise RG calculations, would be of high interest. 

GS acknowledges the financial support provided 
through the European Community's Human Poten- 
tial Program under contract HPRN-CT-2002-00307, 
DYGLAGEMEM. 



where F^{w,u) is a universal |l8l | scaling function: 

F^{w,u)^u{Fp{w)~Fp{w^)) (21) 

where F^{w) — — J_^dw'Fji{w') and je is the Euler 
constant. The exponent A in Eq. ((HJ is defined such that 
lim„_+oo F^{w, u) = 00 (w) with, in the limit T ^ 0: 

47r2 T 

Fcoo{^) = [1 + —^)Fn{w) (22) 

Given the scaling form obtained, and the discussion be- 
low Eq. jnj, it is thus consistent to compare our results 
for the equilibrium exponent zco to its value ob- 
tained in the numerical simulation of Ref. jl4l |. This 
comparison is shown on Fig. 2. We compare the nu- 
merical results both to the low T expansion of zco , 
extrapolated to all temperatures (FRG2), and to the full 
expression of zqq where the value of — rA"*(0) is ob- 
tained from the numerical solution, at finite T - although 
only valid at low T - of Eq. (HHJ (FRGl). Both FRG es- 
timates suggest a rather good agreement, at low T, with 
numerics. 

Finally we can evaluate the non trivial Fluctuation Dis- 
sipation Ratio (FDR) X^^, defined from the connected 
correlation (i.e. such that X^^, = 1 at equilibrium) and 
found to take the form: 

{X^,,)-' = d,,Cl,/iTnl) = Fxiq^t - t'),t/t') (23) 

In the limit t/t' ;» 1, keeping q^{t — t') fixed, one obtains, 
in the low T limit, using (fHll^ : 
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